function Y = F(U)
% F - x方向通量
%   
[rho,u,v,p] = U2Variables(U);
E = U(4);
Y=[rho*u;
   rho*u^2+p;
   rho*u*v;
   u*(E+p)];